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ABSTRACT 

The basic parameters and boundary conditions for 
solution of the equilibrium equations for a viewport in 
a deep submersible are defined.. A likely solution method 
is chosen and investigated in order to determine its 
applicability and to learn more of the behavior and 
inter-relationships of the governing equations. The 
difficulties encountered in the chosen solution method 
are explained and another approach is suggested. 
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INTRODUCTION 



Since Piccard's Trieste,^ - almost all deep 
submersibles have utilized small plexiglass viewports 
for observation and control. The governing structural 
mechanics equations for the viewports have not been 
solved and, therefore, the stress state and predicted 
deflections have not been known. The interaction with 
the pressure hull and its reinforcement has also mainly 
been known from finite element analyses. 

The principal investigations of the behavior of 
plexiglass viewports were experimental tests reported by 
Stachiw. ^ ^ ' 4 ' ^ These tests included short term loading 
of conical, flat, and spherical windows as well as long 
term creep of conical windows. The data from these tests 
provides sufficient information for safe, large safety 
factor design. It does not, however, provide insight 
into the stresses in the window support or detailed 
information of expected viewport performance and useful 
life. 

Since plexiglass is inexpensive, the viewports have 
been greatly over-designed and the present state of the 
art has been adequate. However, for future, more 
sophisticated vehicles with larger viewports, a better 
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understanding of the stresses and deflections in the 
viewports and the resultant interaction with the pressure 
hull will permit a safer and less expensive design and 
greater insight into critical design aspects. 



BACKGROUND 



Most viewports presently utilized in undersea vehicles 
are truncated cones mounted in heavily reinforced openings 
in the pressure hull. They are held in position by 
mounting rings but the primary seal is simply from surface 
to surface contact. The contact surfaces are usually 
polished to ensure a good seal and a 90 degree cone angle 
is utilized. This angle was determined by Piccard to be 
the best design and has been utilized successfully by 
designers since then. The tests by Stachiw have supported 
the opinion that the 90° cone angle is best for acrylic 
viewports . 

It can be seen that the viewport is similar to a 

thick, flat, circular plate with slanted ends. The stress 

distribution is axi— symmetric and therefore circular 

coordinates should be utilized in forming the governing 

equations. For the coordinate system in Figure 1, the 

deformation is symmetrical around the z axis and therefore 

all derivatives with respect to 0 are zero and shear 

stresses x „ and x„ are zero, 
r w 0 z 

Timoshenko and Goodier^ show that the equations of 
equilibrium therefore take the form: 
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Figure 1 

/ 
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( 1 ) 




At present most submersibles have moderate depth 
capabilities and therefore do not need excessively thick 
viewports. The desire for good visibility therefore results 
in a design with a thickness to diameter ratio of usually 



only tested viewports with thickness to inner diameter 
ratios of .631 or less illustrates that this is presently 
the major area of concern. Consequently, the compatibility 
equations can therefore be simplifed by assuming plane 
stress. Although this is not completely accurate, it 
should allow an accurate prediction of the behavior of 
almost all desirable viewport designs. 

g 

The compatibility relations therefore take the form: 



less than .5 for the viewport. The fact that Stachiw 
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The strain components, therefore, become: 
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Equations (1) and (2) can be satisfied by expressing 

g 

the stresses in terms of a stress function </> . The stresses, 
therefore, become: 
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subject to the condition that: 



2 2 

V V 4> = 0 



( 3 ) 
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These equations apply to a viewport in an undersea 
structure or pressure facility provided the appropriate 
boundary conditions are met. The first and most obvious 
boundary condition is that the normal stress on the 
pressurized face is equal to and opposing the external 
pressure. Obviously, the condition is reversed on the 
unpressurized face where the normal stress is zero. The 
shear stress on the flat faces is also zero since they 
are free. 

The boundary conditions on the slanted ends are less 
obvious. In reality there will be a contraction and 
rotation of the mounting reinforcement with pressure. 
However, since the reinforcement is very thick steel 
with a high modulus of elasticity, it can be safely 
considered as being rigid compared to the plexiglass. 
Since the viewport is not physically attached to its 
mounting, it is free to slide along the mounting opposed 
only by friction and the shape of the mounting. 

Two boundaries on the behavior of the viewport at 
its mounting can be specified. The viewport can be 
considered to be rigidly fixed to its mounting or it can 
be allowed to slide freely. In the first case, it seems 
that the stresses produced will be higher than the actual 
situation and the displacements will be lower. In the 
second case, it seems that the stresses produced will be 
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slightly lower than the actual situation and the displace- 
ments will be higher. 

The first case is satisfied if all displacements 
are zero at the boundary, i.e. 

D n H 

u=w=0atr= — - (j + z) tan 6 



The second case is satisfied if the displacement normal 
to the mounting is zero at the boundary, i.e. 

D fl H 

u cos 0 + w sin 0 = 0 at r = — - (— + z) tan 0 

2 2 



In seeking solutions of equations (1), (2), and (3), 

Timoshenko and Goodier considered a polynomial form of 
the stress function <j> which also satisfied the equation: 




+ 



2 6 $ 
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R 2 <5 Y 2 
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This resulted in particular solutions of the form: 



♦n = a 0 lz 



n 



n(n-l) „2 „(n-2) ,n(n-l) (n-2) (n-3)„4 „(n-4) 
2 (2n-l) z 2-4 (2n-l) (2n-3) 



and additional solutions of equations (1) and (2) found 

2 

by multiplying the particular solution by R where 

2 2 2 

R = r + z . Consequently, the following possible 



solutions can be found: 



% = a o 



4>, = 



a l z 



* 2 = a 2 [2z 2 - r 2 ] 



<J> = 



a 3 [ 2z 3 - 3r 2 z] 



<f> 4 = a 4 [8z 4 - 24r 2 z 2 + 3r 4 ] 



<l> = 



k 2 3 4 

a c [ 8z D - 40r + 15r z] 



<j> = a [^fz 6 - 40z 4 r 2 + 30z 2 r 4 

6 6 3 



5 6, 

3 r ] 



and 



*2 


= b 2 [r 2 + 
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Timoshenko and Goodier considered a flat circular 
plate loaded with a uniform pressure on one face. The 
boundary conditions were: 
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x = 0 at z = 
rz 2 



a = 0 at z = 

Z 2 



° z - - p at z = -f 



/ 



These are identical to all but the last boundary condition 
for a viewport. 

By assuming a solution form of : 



* ,16 6 4 2 , 2 4 56 , 

<P = a, (—~—z ~ 40z r + 30z r - -^-r ) 

D J j 



I L-. / rt 6 T ^ 4 2 ri -l 2 4,^6l 

+ b, (8z - 16z r - 21z r + 3r ) 

6 



+ a 4 (8z 4 - 24r 2 z 2 + 3r 4 ) 



+ a^(2z 2 - 3r 2 z) 



+ b^ (z 2 + zr 2 ) 



they found that 



a 6 = ( ~ — r^ } P 

0 3520H J 



b fi = ! —*) ? 

352H 
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a 4 ( 384H } P 



a = ( ~ 1+5 v ) p 

3 '72+60 v r 



b 3 ( 24+20 v 



) P 



Therefore the stresses become: 



a = P[(2 +v)4 - 

r H H 



— — ] 
4 H J 



a = P[- 



2H 



: + 2 £ _ I] 

3 + 2 H 2 J 



T rz " 3pr <V " z2) 



Using the equations for strains given previously we 
find that: 

1+v 5 2 <p 

u E Sr6z 

w = [ 2 ( 1- v ) — — ^ + 2 ( 1— V) i || + (l-2v) -M^] 

Jbj r z r or j z 

or Sz 

Consequently , 
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z r+ ( — _ ) zr + (-t-=j) zr+ ( 
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1+v r ,-l-V x 4 , , 3 v x 2 2 , ,3-3v x 4 , , 3 x „2 , , 3 x _2 

w = ” 1 ( l? ,z +< ^ )z r +< I^ )r + < ® )z +( ' ™ ,r 



+ ( 12+10v )zJ P 



The final boundary condition for the first case for 

D q jj 

the viewport requires that u=w= 0 atr= — - (j + z)tan 0 
The flat plate solution, however, gives 



1+v ,-2+v. 4 . fl ^ /3-3V. 4 3 a 

u = - 7 =— P [ { — ^—) z tan 9 + ( ~) z tan 0 

E H 3 4H 3 



,2-v x . 3 . ,-9+9 v. A 3 . 2 a 
+ ( o) A z + ( ■=— ) A z tan 0 

2H 8H 



9-9%, 22 3 2 

+ ( ^ ) A z tan 6 + (-j-rj) z tan 0 

16H 3 4H 



+ < 32H^ ) + ^ AZ + 2 tan 6 



+ ( 24+20v ) A ' 0 



_ 1+V , ,-1-v. 4 , - 3v. 4 . 2, 

W = — pr- p [ ( — 5—) z + ( — =-) 7, tan 0 

2H 211 



+ ( 3 3 ^ ) z 4 tan 4 0 + (— ^-) A 7 ^ tan 0 
16ir 2H J 
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+ <=2±|v, Az 3 tan 3 e + , 3 , 2 2 + ,jv J a 2 2 2 
8H 3 8H 3 



+ (■^•) z 2 tan 2 0 + { - A 2 z 2 tan 2 0 

8H 32H 3 



+ ^12+lOv 5 2 + { 8H ) Az tan 9 + ( 32r 3 } A 2 tan 9 



+ (~t) A 2 + ) A 4 

32H 256H 3 



The final boundary condition for the second case for 

the viewport requires that u + w tan 0 = 0 at 
D 0 H 

r = [ — - (— + z) tan 0 ] . The flat plate solution, 
however, gives 



, n 1+v _ r 5+v. a 4 , 3+3v. 3,-. 4 

u + w tan 0 = — P [ ( — ~— ) tan 0 z + ( ^-) tan 0 z 

E 2H 2 4H 3 



,3-3v. 5 q 4 .2 — v. . 3 , , — 9 — 3v. . 2„ 3 

+ ( =• ) tan 0 z + ( — 5 -) Az + ( —) Atan 0 z 

16H 3 2H 3 8H 3 



,-3+3v> 4 q 3 , ,9-3v. .2, Q 2 , , 3, . 3 Q 2 

+ ( 5 —) Atan 0 z + ( A tan 0 z + (-• ott) tan 0 z 

8H 3 16H 3 8H 
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+ (A) tan0 z 2 + ( 




+ (- l)Az + (^r)Atan 2 9z + ( 



12+10 



-5 



) tan 0 z 



+ ( 3 ~ 3 . ~y .) A 4 tan 0 + (- 
256H -3 




±1 ) A ] 

24+20v' J 



where A = (Dq - Htan 0 ) . 

Consequently, it can be seen that if these boundary 
conditions could also be satisfied, solutions would result 
which would be dependent on the material properties v and 
E as well as the dimensions D^, H, and 0 and the pressure 
P. Since plexiglass is a viscoelastic material the 
properties v and E will be slightly non-linear and will 
be functions of time. In practice, all viewports are 
designed to remain well within the elastic range since 
vision would be greatly distorted even before reaching 
the non-linear response. Consequently, a solution 
assuming that the material is linearly elastic will be 
adequate since the variation in properties of different 
batches of plexiglass makes more exact calculations of 



little use. 



SOLUTION WITH SIXTH ORDER 



POLYNOMIAL ASSUMPTION 



In order to satisfy either of these additional 

boundary conditions it is apparent that supplemental sixth 

order terms will be required. Since no other terms are 

possible using the method of Timoshenko and Goodier, 

other forms must be found. The desired forms must 

satisfy equation (3) and two possibilities therefore 

2 2 2 

become apparent. Either VV(p = 0 or Vtp = 0 will 
satisfy equations (1) and (2) . 

Frobable desired terms include sixth, fourth, and 
third order terms. General terms can be constructed as 
follows : 



<j> = az^ + 
<J> = ez + 
<P - pz ^ + 
In order to satisfy 

a = 



,42 
bz r + 



-22 
f z r + 



2 4,6 
cz r + dr 



4 

gr 



2 

qzr 

V^V^cp = 0 we find: 

12b + 8c 
45 



3b + 8c 
72 
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Additional terms can be constructed as follows: 
<f> = hz^ + iz^r^ + jz^r^ + kr^ 

<p = lz^ + mz^r^ + nr^ 

2 2 
<f> = sz^ + tr 

2 

In order to satisfy V <f> 

h = 
i = 
j = 



= 0 we find: 




24k 

-18k 



m = -8n 




The boundary conditions for the combination of the 
flat plate solution and the additional terms therefore 
become : 



T 



a 



rz 



z 



0 at z = ± ^ 

n . H 

0 at z = 2 
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0 4- H 

a = -P at z = - 7T 

z 2 



u = w = 0oru + w tan 0 = 0 at r '= ^-[D^ 



(H+2 z) tanQ] 



The boundary conditions for the additional terms alone 
therefore become: 



T 

rz 



0 at z 




a = 0 at z = ± ^ 
z 2 



u + U,, , , , = w+w £1 , n 

flat plate flat plate 



0 at r= 




z 



or 

u + wtan0 + 



D 0 

(u+wtan 0),, . , . =0 at r= — 

flat plane z 



,H 



(^- + z) tan0 



+ z)tan0 



For the additional terms to satisfy the first 
boundary condition we find: 

T rz = 0 = (8c - 8v c - 6vb - 144k) H 2 r 

+ (-6b + 6vb - 16c + 8vc + 144k) r^ 



+ (32g - 32vg - 4vf + 32n) r 
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Since the constants b, c, k, g, f, and n are independent 
3 

of r, the r and r terms must each equal zero. Therefore 
144k = 6b - 6vb + 16c - 8vc 

and 

32n = (6b + 8c) H 2 + (-32g + 32vg + 4vf) 

Considering the second boundary condition we require 

that: 

a = 0 = (-4b - c) (±H 3 ) + ( 8q - 4vq + 6p - 6vp) 

Z j 

Therefore 



6p (1-v) = -8q + 4vq 



and 



b = - 



4 

3 C 



D 



At the boundary r = - 2 - - 



(^- + z) tan 6 



u = — [ (-16b + 8vb - c + vc) z^ tan 0 



+ (6b - 6vb + 8c - 8 v c) z 4 tan 3 0 



+ (8b - 4vb + c - v c) Az 3 
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+ (-9b + 9vb -12c + 12vc) A 2 z 2 tan 2 0 + (6b + 8c)H 2 z 2 tan0 
+ (|b-|vb+ 6c - 6vc)A 2 z 2 tan6 

+ (-32g + 32vg - 4f + 4vf)z 2 tan6 

+ (-|b+ |vb- c + vc) A 3 z + (-3b - 4c)H 2 Az 

+ (16g - 16vg + 2f - 2vf)Az + (-2q) z tan 6 + (q) A ] 
and 

1+V r /ylk _i_ A W . 16 ,16 . 4 

w = — ^ — [(4b + 4vb + — c + — VC) z 

Hi 6 6 

+ (-12vb - 16v c) z 4 tan 2 0 + (-|b+ |vb - 2c + 2vc) z 4 tan 4 0 
+( 12 vb + 16vc)Az 2 tan0 + (3b - 3 vb + 4c - 4vc)Az 2 tan 2 0 

+ (-3V b - 4 v c) A 2 z 2 + (-|b+ |vb- 3c + 3vc) A 2 z 2 tan 2 0 
+ (-6b - 8 c ) H 2 z 2 + (4vf + 32vg)z 2 + (3b + 4c)H 2 z 2 tan 2 8 



+ (16g - 16vg + 2f - 2vf ) z 2 tan 2 0 + (-3b -4c)H 2 Az tan 0 
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+ (jb - I vb' + c - ^c)A^z tan 0 
4 4 

+ (~16g + 16vg - 2f + 2vf)Az tan 9 + (6p - 12vp + 8q-8vq) z 

+ { ~j2 h+ 72 Vb " | c+ 7 vc ) a4 + (fb+ c)H 2 A 2 

+ (4g - 4vc + if - i vf) A 2 ] 

Substituting for p and b we find: 

u = [ (-32g + 32vg - 4f + 4vf)z 2 tan 6 

+ (16g - 16vq + 2f - 2vf) Az (4) 

+ (-2q) z tan 0 + q A ] 

and 

1+V 2 22 

w = — - [ (4vf + 32vg) z + (16g - 16vg+2f - 2vf)z tan 0 

4v 

+ (-16g + 16vg - 2f + 2vf)Az tan 0 + 

+ ( 4g - 4vg + if- i vf) A 2 ] 



z 



(5) 
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Consequently, it is apparent that neither of the 
final boundary conditions can possibly be satisfied for 
the third or fourth order terms using this solution form. 
This occurs because the first boundary condition causes 
the fourth order constants, f and g, to cancel in the 
second boundary condition. This produces a relationship 
between b and c which prevents them from contributing 
to the displacements u and w. However, the addition of 
higher order terms may allow all boundary conditions to 



be met. 



SOLUTION WITH INFINITE ORDER 



POLYNOMIAL ASSUMPTION 

In order to investigate this possibility we consider 
a stress function of order n. Therefore 



, n , , n-2 2 , n-4 4 

<f> = a z + bz r + cz r + .... 

n n n n 



The terms contain only even powers of r since only even 

2 2 2 

powers will satisfy equation (3). Therefore V <f> and V V <J> 
may be expressed as: 

V 2 <f> = [2-2b + n(n-l)a ] z n 2 + [4*4c + (n-2) (n-3)b ] z n ^r 2 

r n n n n 



+ [6*6d + (n-4)(n-5)c ] z n ^ + .... 

n n 



and 

V 2 V 2 4> = [2*2*4. 4c +2*2*2(n-2) (n-3)b 
r L n n 

+ n(n-l) (n-2) (n-3) a n ] z n_4 + [ 4 • 4 • 6 • 6d n +2 • 4 • 4 (n-4) (n-5)c n 

+ (n-2) (n-3) (n-4) (n-5)b ] z n-6 r 2 + [ 6 • 6 • 8 • 8e +2 • 6 • 6 (n-6 ) (n- 7) d 

n n ' n 

+ (n-4) (n-5) (n-6) (n-7) c ] z 11 8 r 4 + .... 
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2 2 

Consequently, in order to satisfy V V <p = 0 we find: 



= _ 2 (n-2 ) (n-3) . _ n(n-l) (n-2) (n-3) 

n 4.4 n 2 • 2 • 4 • 4 n 



_ 3(n-2) (n-3) (n-4) (n-5) w . 2n (n-1) (n-2) (n-3) (n-4) (n-5) „ 

n 4 • 4 • 6 ♦ 6 n 2. 2. 4. 4. 6. 6 n 



_ _ 4 (n-2) (n-3) (n-4) (n-5) (n-6) (n-7 ) , 
n ' 4 • 4 • 6 • 6 • 8 • 8 D n 



3n (n-1) (n-2) (n-3) (n-4) (n-5) (n-6) (n-7) 

2*2*4*4*6*6*8*8 a n 



etc. 

. 2 
Using Greek letters to distinguish terms, V <J> = 0 is 

satisfied if: 



B = - 



n 



n (n-1) 
2-2 



a 



n 



y _ n (n-1) (n-2) (n-3) 
y n 2-2 - 4 - 4 n 



^ _ -n (n-1) (n-2) (n-3) (n-4) (n-5) 

n 2- 2- 4- 4- 6- 6 n 



etc. 

The shear stress at any point can therefore be 



represented as: 
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T rz = [(-4 + 2v) { (n-2) (n-3)b n ) + (-1+v) - - -y- 2 - - ~ >a n 

+ n (n~l) (n- 2 ) (n-3) a ^ z n ~ 4 r + [ (3 y) { ilLzil ( n ~ 3 ) (n-4) (n-5) 

2 n 4 

J. /o ^ j-n(n-l) (n-2) (n-3) (n-4) (n-5) 

+ (2-v) t JTJ71 }a n 

n(n-l) (n-2) (n-3) (n-4) (n-5) .. , _n-6„3 

— — — CX Z IT 

2*2*4 n J 

+ [(-4 + v) { (n ~ 2) (n ~ 3) (n ~ 4) (n ~ 5) (n ~ 6) (n ~ 7) ] b 

4*4*6 n 

+ (_ 3 + V ) [ n (n-1) (n-2) (n-3) (n-4) (n-5) (n- 6 ) (n-7) 

2 • 2 • 4 • 4 • 6 n 

, n(n-l) (n-2) (n-3) (n-4) (n-5) (n- 6 ) (n-7) , „n-8„5 

2 • 2 • 4 • 4 • 6 n JZ r 

etc . 

The axial stress can also be represented at any 
point as : 

a = [ ( 2 -v) 2 - 2 (n- 2 )b + ( 1 -v) n(n-l) (n- 2 ) a 

^ n n 

- n(n-l) (n- 2 ) a n J z 11-3 

+ t (~3+v) (n-2) (n-3) (n-4)b 4- (-2+y) ^-UirjI(PJl3)_( n- 4X ■ 

n 2-2 n 

n(n-l) (n-2) (n-3) (n-4) , „n-5„2 

+ 2^2 V* r 



